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1. Introduction 

The asymmetric simple exclusion process (ASEP) is one of the most studied models 
far from equilibrium It has been used to describe various problems in many 

fields of research, such as biopolymerization and traffic flow. A lot of analytical results 
exist, both for open and periodic boundary conditions. The model describes a particle 
system on a chain with hard core exclusion. Particles are allowed to hop one site 
to their right, supposed that it is empty. The usual dynamics is in continuous time: 
the random-sequential update [Tj. Moreover, discrete-time update schemes have been 
studied: backward- and forward-ordered sequential updates, site-oriented jl] as well 
as particle-oriented [5], sublattice-parallel updates j6j [7] and fully-parallel dynamics 
[HI HE UHj (for an overview, see [TT]). We analyze another update scheme, that has 
originally been introduced in a two-dimensional cellular automaton describing pedestrian 
dynamics O HB] . 

A typical situation encountered in pedestrian dynamics is the motion along a 
corridor. This is modeled as a strip of finite width W in ^/-direction of a cartesian 
coordinate plane and length L in x-direction divided into square cells. Each cell can be 
in one of two possible states, i.e. either occupied by one of the N pedestrians, or empty 
(hard core exclusion rule). Assuming periodic boundary conditions in x-direction and 
impenetrable walls in ^-direction implies that the density of pedestrians p = N/ (W-L) is 
constant. A pedestrian can move to neighbor cells with different transition probabilities, 
depending on the direction and the particular neighborhood. The preferred direction is 
along the x-axis. A parallel update, in which all particles are updated simultaneously, 
would lead to conflicts, in which more than one person tries to access the same cell |15| . 
Hence, another way of updating was chosen, the so-called shuffled dynamics: At each 
(discrete) timestep, the order in which particles are allowed to move is determined by a 
random permutation. Since the theoretical implications of this update procedure have 
not been considered before, we here investigate the one-dimensional limit first. In the 
language of pedestrian dynamics this corresponds to a very narrow corridor with W = 1 
such that 'side-by-side' motion or overtaking is not possible. The model then becomes 
equivalent to the ASEP with shuffled update. 

The difference between this update scheme and the random-sequential update 
mentioned above is analogous to an urn problem with and without replacement: If 
one imagines an urn with N balls, numbered 1, 2, . . . , N, the random-sequential update 
can be realized by choosing a ball at random, updating the particle with the ball-number 
and replacing the ball into the urn afterwards. In the shuffled update, one chooses a ball 
and updates the corresponding particle without replacing the ball. Then, one chooses 
the next ball, and so on, until the last particle is updated. Then the urn is refilled and 
the procedure is repeated in the next timestep. 

After giving a precise definition of the model, approximative formulas for the steady 
state distribution and fundamental diagram are derived. By mapping onto generalized 
zero-range processes it is shown that the ASEP with shuffled update does not have 
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a product measure steady state for general parameter choice. Some applications and 
generalizations are given afterwards and concluding the results are discussed and some 
interesting directions for further research are given. 

2. ASEP with shuffled dynamics 

In the following we give a more formal definition of the ASEP with shuffled dynamics. 
Consider a one-dimensional lattice with L sites and periodic boundary conditions. Each 
site may either be occupied by one of the N particles, labelled % = 1,2, ...,N, or it 
may be empty. Therefore the particles are distinguishable. In each discrete timestep 
a random permutation 7r(l, ...,N) of the particle labels equals the update sequence. If 
the right neighboring cell is empty, the relevant particle moves one site to the right with 
probability p, if it is occupied, the particle stays in its cell. 

Figure Q shows a part of a large system consisting of six cells and four particles 
(numbered 1,2,3,4 without loss of generality), at time t (left) and t + 1 (right). The 
drawn update sequence is . . . , 3, . . . , 4, . . . , 1, . . . , 2, . . ., where the ellipsis indicate that 
other particle numbers belonging to different clusters (units of neighboring particles) can 
be chosen between (for the cluster depicted, only the relative positions of the numbers 
of its particles in the sequence are of interest). Particle 3, chosen first, can not move 
since the cell in front is occupied by 2, and similar for particle 4. Considering the case 
p — 1, particle 1 then moves deterministically to the right. Then 2 also moves, because 
it was drawn after 1. Although particle 4 is drawn after 3, it can not move, since both 
were drawn before 1 and 2. 



4 3 2 1 4 3 2 1 




Figure 1. Shuffled update of a cluster consisting of four particles numbered from right 
to left. The drawn sequence is 3, 4, 1, 2 and p = 1. 

As pointed out in the Introduction, the shuffled update is different from the random- 
sequential dynamics which is generically used for the ASEP. Whereas the latter describes 
stochastic processes in continuous time, the shuffled dynamics combines elements of 
discrete updates and dynamics in continuous time. It is discrete in the sense that there 
is a well-defined timestep during which each particle is updated exactly once. On the 
other hand, the order of updating the particles is not fixed. E.g. it may happen that a 
specific particle is updated last during a timestep and first during the next one! This 
can actually be a problem for applications because it is difficult to identify the timestep 
as a kind of reaction time (as it is natural for the case of parallel dynamics jH|). 

Despite this important difference, the two updates share certain similarities. In 
contrast to the ordered updates with fixed order, they do not have a deterministic 
limit, even for hopping probabilities p — 1. However, in the random-sequential case the 
dynamics depends on p only in a trivial way, since by rescaling time always p = 1 can 
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be chosen. This is not possible in the shuffled case. We therefore expect a non-trivial 
p-dependence of the results, as in the other discrete-time updates. 

3. Steady State Distribution 

Using a particle-oriented representation |H] , the state of a system with N particles 
at time t can uniquely be specified by \ni(tY ni ^\ . . . ,riN(tp nN ^}. Here rij represents 
the state of the j-th particle (i.e. the number of empty sites in front) and ttj its update 
number (the position of the j-th particle in the ordered update sequence). At each 
timestep, the permutation operator V(t) generates the new random sequence as follows: 
V(t)\ . . . , n^i)^*^, . . . , ) = | . . . , rii{t)^ t+V) \ . . .), i.e. without changing the variables 
Wj(t). The complete update of the system t — > t + 1 can be described by the action of 
the transfer- matrix hP^hS 2 ^ . . . hS N \ being an ordered product of the local operators hP* 
acting on the state of the particles <Jj—\ and Gf h^ = h^_ la . = l—p{a (T:j a^ IJ] —a\._ 1 a (Tj ) 1 
where <x,- is the particle with update number j. + The master equation reads 

(P\V(t) • . . . h® ... | ... , m{t)^\ • • •) = (P\ • • • , n t (t + i)(^ +1 )), . . .), (1) 

where (P\ . . . , rij(£), . . .) denotes the probability of the configuration | . . . , rii(t), . . .}. In 
the steady state, equation (JTJ) simplifies, because all time-dependencies vanish. Since the 
permutation operator generates random sequences (each of them with probability 1/iV!) 
one can calculate the probability with which a particle moves in a certain configuration. 
The n-th particle of a cluster has moved at time t + 1, if and only if at least the first n 
particles are chosen in the order from the right to the left (with probability 1/n!) and 
if they all moved (with probability p n ). The probability that exactly / particles of a 
cluster of length n (with n > I) move is given as the probability for / particles to move, 
as calculated above, minus p k+1 /(k + 1)!, the probability for the (/ + l)-th particle to 
move. We arrive at ui(n), the probability for I particles leaving a cluster of length n 
(where < I < n): 

(and otherwise). The Heaviside step function ensures that in the case / = n the right 
term vanishes. Note that although the particles in the ASEP are updated in N steps, 
it can be considered as if all clusters were updated in parallel. 

We now leave the operator notation and write the steady state probability simply 
as P(ni,n 2 , . . . ,n^). We assume that this probability can be written as a product of 
the probabilities P Ui to find particle % with rii holes in front, i.e. 

P(m, n 2 ,..., n N ) = P ni P n2 ■ ■ ■ Pn N - (3) 

This constitutes the so-called car-oriented mean-field (COMF) theory, successfully 
applied previously to traffic flow models [22123 • Note that equation © implies that only 

+ For a detailed explanation of the used formalism, see 0- 
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correlations between neighboring particles are taken into account. In the thermodynamic 
limit, the master equation for the steady state reads: 

p =p -(p- g) (l - Po)+pgPu (4) 
Pi = (p- g) (i - Po) + (pg + Pg) Pi + pgP%, (5) 
P n = pgPn^i + (pg + pg) P n + pgPn+i, Vn > 1. (6) 

Here Pk (k = 0, 1, . . .) is the probability for an arbitrary particle to have k holes (a 
hole-cluster of length k) in front and g = 1 — g is the probability that the particle in 
front moves. As an example we explain the equation for Pi in the following: A gap of 
1 hole can arise from three different processes corresponding to the three terms on the 
r.h.s. of (0): If the particle has 2 holes in front (probability P2) the gap decreases by 1 if 
the particle itself moves (probability p) , but the preceeding one not (probability g) . The 
gap remains unchanged (probability Pi) if either both particles move (probability]?^) or 
if both particles do not move (probability pg). If the particle has its preceeding particle 
directly in front (probability P ) the situation is more sophisticated: The particle has 
one hole in front afterwards if and only if the particle in front has moved, but the particle 
itself not. This probability depends on the position of the particle in the cluster to which 
itself and the preceeding particle belong. The probability P{k) for the particle to have 
exactly k — 1 particles directly in front is approximated by 

P(fc)=P fc - 1 (l-P ). (7) 

This probability (JJJ) has to be multiplied with u k _i(n) = rfczyi ~~ fp f° r k < n, 
obtained from equation (|2j). Finally we have to sum over all possible k. This yields 
(1 — Po) Y^k=2 ( (fc-i)! — fr) ^0 w hich can be rewritten as the first term on the r.h.s. 
of (0) as one can check easily if one interpretes g as the hopping probability for a particle, 
averaged over all possible numbers of particles in front, i.e. 

fP, for P = 0, 

l — ^ (exp(pP ) - 1) , for P > 0. (§) 

•M) 



00 u 
P 

g 



fc=i 



Using this, the system of equations (fl]) - can be solved by generating functions 
|2*2"1 12*3] . We obtain an implicit expression for P that has to be solved numerically in 
general, since g and g depend on P via (jHJ): 

Po = tiP-P)-lPP-p)9 (9) 

ppg 

D ip-g)(l-Po) f Pg 



\/n > 0. (10) 

pg \pg; 

Again, the abbreviations p — 1 — p, etc. were used. Since g is the probability that an 
arbitrary particle moves, it equals also the average velocity v of the particles. This is 
related to the flow J and yields the so-called fundamental diagram 

J(p) =vp = gp. (11) 
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In the (partially deterministic) case p — 1, Pq becomes 




(12) 



as one can check easily. Note that p = 1 and p < 1/2 is the only case in which Pq can 
vanish. For p = 1, the probability Pq is completely determined by the fact that the 
first particle of a cluster moves deterministically and all the other particles have smaller 
hopping probabilities due to the shuffling. For densities p < 1/2 this implies that any 
state which consists only of clusters of size 1, i.e. separated particles, is stationary. Hence 
the probability to find a particle directly in front vanishes and we have Pq = (Fig. |2J). 
For densities p > 1/2 clusters are formed that are separated by exactly one hole, i.e. 
there are only isolated empty cells in the steady state. It is easy to verify that no pairing 
of holes can happen, since from the point of view of the holes they jump at least one 
site backwards and at most to the end of the cluster. As a consequence, P n = for 
n > 2. Therefore L — N clusters exist, which is then also the number of particles having 
exactly one hole in front. Thus we obtain Pi = (1 — p)/p and Pq = 1 — Pi = (2p — l)/p 
in this density regime. These results are exact for p = 1 and are reproduced by COMF 
(see ©-(HI). 

The calculation of the fundamental diagram using (|TT|) requires the knowledge of 
g. This in turn depends on the cluster length distribution P(k) which is not known 
exactly. Using the approximation (fTTj) . the flow-density relation is explicitly given by 



The fundamental diagram (|13|) shows a strong asymmetry with respect to p = 1/2. 
For densities p < 1/2 each particle can move independently and deterministically for 
p < 1/2, since every particle has at least one hole in front, exactly as in parallel updating 
|23j . If the density is increased to values greater than 1/2, L — N clusters of nonvanishing 
length are formed from which the rightmost particle can move deterministically and 
since the probability to find such a particle in the system is given by 1 — p, they 
add exactly this value to the flow, as in the usual parallel update. Consequently, this 
result can be obtained from ()13|) by a first order Taylor expansion. The contribution 
of the other particles to the flow depends exponentially on the ratio of these particles 
(p — (1 — p))/p — (2p — l)/p. This yields the curvature in the fundamental diagram in 
the high density regime. In contrast to the parallel and random-sequential dynamics, 
the shuffled update is not particle-hole symmetric (Fig. |2j). 

In this deterministic case clearly a free- flow phase (p < 1/2) and jammed phase (p > 1/2) 
can be distinguished. They are separated by a phase transition at p c = 1/2. The 
probability Po(p) constitutes a kind of order parameter. It vanishes exactly for p < 1/2 
and increases continuously for p > 1/2 (Fig. |2J) indicating a second order transition. 
This is similar to the deterministic limit of parallel dynamics [26] . With decreasing 




forp< 1/2 



forp> 1/2. 



(13) 
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Figure 2. The probability Po(p,p) (left) and the fundamental diagram J(p,p) (right) 
for p — 1 (top) and p — 0.5 (bottom). Depicted are the results from COMF (lines) and 
from computer simulations (squares) for a system consisting of L = 500 cells. Note 
that Po(p,p = 1) is exact. 



hopping probability the asymmetry of the fundamental diagram becomes smaller. In 
the limit p — > we obtain the result for random-sequential update, J = pp(l — p). For 
p < 1 the free-flow and jammed regimes are no longer separated by a phase transition 
as can be seen from the smooth behaviour of -Po(p) (Fig. |2j) . 

We just mention [21] that the assumption of a factorized steady state in a finite 
system yields less agreement, but the analytical results converge fast to the result of 
the thermodynamic limit, derived here. In case p = 1, for example, the largest possible 



cluster has length 2N — L + 1 and the term 



2p-l 



exp 



2,0-1 ' 
P ' 



- 1 



in (fT3j) has to be 



replaced by the generalized hypergeometric function 2 .F2(1,L — 2N; 2, 1 — N; 1) [25] . 
which overestimates the results for small L, indicating the mean-field result is not exact, 
at least for finite systems. 



4. Mapping onto a generalized zero-range process with parallel dynamics 

The ASEP can be mapped onto a model with multi-occupation of sites, the generalized 
zero-range process (GZRP). Unoccupied cells, i.e. holes, of the ASEP become the sites 
of the GZRP and the particles between the (i — l)-th and the i-th hole of the ASEP 
(the particles of the (i — l)-th cluster) now all occupy the i-th site of the GZRP and are 
referred to as mass mi (see Fig. Ola) and b)). We have L — N sites and masses and still 
N particles. At each timestep, particles are chipped off the mass trii located at site i 
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(i = 1, . . . , L — N) with transition probabilities ui^rrij) and are moved to site i + In 
the ASEP with shuffled dynamics the probabilities ui(m) are given by (J2J). Note that 
although the particles in the ASEP are updated in N steps, it can be considered as if all 
clusters were updated in parallel. Thus indeed the ASEP with shuffled update is mapped 
onto a GZRP with parallel dynamics. In the case p — 1 and N < L — N, we have seen 
that in the ASEP the probability to find two neighboring particles vanishes. For the 
GZRP this implies that each cell is either occupied by exactly one particle, or it is empty. 
Thus the GZRP is in fact a usual ZRP in this case and its steady state distribution can be 
written as a product measure. The authors of j^j have derived a necessary and sufficient 
condition for the existence of a factorized steady state distribution in a broad class of 
one-dimensional mass transport models, including the GZRP. This is possible if the 
transition probabilities Ui(n) can be written as a product viw n -ij Y^i=o v i w n-h where v\ 
and w n -i are functions that depend only on I and n — l, respectively. They also derived a 
more direct test in . For continuous-time dynamics a related condition has been found 
recently for processes in arbitrary dimensions [201 - T ne transition probabilities (J2J) of 
the corresponding GZRP do not satisfy this condition, for general p and p and therefore 
the steady state does not factorize, i.e. P(mi,m 2 , . . . ,m L _ N ) 

the ASEP, this implies that the particle-cluster probabilities do not factorize. However, 
a factorization into hole-cluster probabilities as assumed in COMF is not excluded. 

5. Mapping onto a zero-range process with shuffled dynamics 

While in the GZRP an arbitrary fraction of particles is allowed to chip off a mass, in the 
usual zero-range process (ZRP) [THj at most one particle (say the topmost) can leave a 
mass during one timestep. The ASEP can be mapped onto a ZRP (T7JHH] by identifying 
the particles of the exclusion model with the sites of the zero-range process, see Fig. El 
a) and c). The holes between the i-th and the (i + l)-th particle of the ASEP form the 




Figure 3. Mapping of the a) ASEP onto b) GZRP and c) ZRP for L = 10 and N = 6. 
The arrows indicate a possible local transition with probability iti(3). 

i-th mass of the ZRP. Thus the ZRP has N sites and L — N particles. Note that the 
particles in the ZRP hop to the left, as the holes in the ASEP. While in the ASEP the 
particles were shuffled, in its equivalent ZRP the sites are shuffled. Zero-range processes 
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are known to have a factorized steady state distribution fT7\ ITS] . This result holds for 
time-continuous dynamics and for parallel and ordered-sequential updating 5J as well. 
For the shuffled update no general results are available yet. The COMF approach in 
Sec. 3 (eq. ©) describes a factorization in the ZRP-picture. For p = 1, the result (fT2"]) 
is exact. However, it is not clear whether the fundamental diagram is exact or not, since 
it depends also on g which we do not know exactly. 

It can easily be seen that this factorization does not become exact for general p, 
by considering the master equation for the probability P(0 Af_1 ,M), i.e. one cell of the 
ZRP having occupation number M := L — N , followed by iV — 1 empty sites (a cluster 
of N particles followed by L — N holes in the ASEP) : 

(N\ N—l y. 

'' fc=i ' 

with p — 1 — p. Inserting the product measure ansatz (jSJ) and performing the limit 
iV — ► oo gives 

| = ^(e'-D. (15) 

Now, inserting the expression for P n , n > 1 from eq. (|T0|) and using (JHJ) leads to the 
constraint 

°= (cP - 1) (^- i ir)- p ' (16) 

This condition is only fulfilled in the limit P — 1 (i.e. for p = 1) and also for 
p = 0. Denoting the r.h.s. of (fT6|) by £(p, -Po) one can see that £ increases almost 
linearly for decreasing P - 111 the limit Po ~ > it becomes maximal, taking the value 
0) = ((2 — p)e p — 2 — p)/2p, which is, for p^O distinctly different from zero. Thus 
one can conclude that the COMF theory gives a good approximation but not the exact 
result, for general p. In other words, the ZRP with shuffled update does not have a 
product measure steady state. 



6. Truncated processes 

In Sec. El we have seen through mapping onto a GZRP that the shuffled update can also 
be interpreted as a cluster dynamics. The probabilities that I particles leave a cluster 
of length n are given by (J2J). These probabilities decrease rapidly with increasing 
Therefore it is natural to consider as an approximation to the original dynamics models 
in which, irrespective of the cluster length, at most Z max particles are allowed to leave a 
cluster at any given timestep. This truncation leads to a GZRP defined by probabilities 

^ ax H^) = ^-^yy^ma X -0, (17) 

with < I < min(7 max ,n) (and otherwise). The case / max = 1 corresponds to the 
standard ZRP with hopping probability u(n) = p, i.e. the ASEP with parallel dynamics. 
Thus, the fundamental diagram is given by J\(p,p) — | f 1 — a/1 — App(l — p)) |E]- 
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In the case Z max 
we have P(m l5 m 2 ,. 
obtain 



2 the factorization condition of ^Sl E] is satisfied for p = 1 and 

For the fundamental diagram we 



mL-N 



P P 

1 mi 1 m,2 



P 



m L -N- 



MP>P = !) 



1 



p+ v /p2_ 2 (l-p)(2p-l) 



for p < 1/2, 
for p > 1/2. 



Fig. shows the results for Z max = 1, 2 in comparison with the shuffled update 
l max = oo. The steady state probabilities for l max = 2 depend on the number of 
particle-clusters with length 1. This is not correctly reproduced by COMF theory which 
assumes a factorization into hole-cluster probabilities Q and yields that all stationary 
configurations for fixed particle number N are equally probable for p — 1 and p > 1/2: 
P(ni, n 2 , • • • , n N ) = p^-LpL-N Therefore 

one can conclude that the steady state 
distribution for the truncated model with l max = 2 and p = 1 factorizes into particle- 
cluster probabilities but not into hole-cluster probabilities. 

For Z max = 3,4,... the curves converge quickly to the shuffled curve, but they no 
longer fulfill the factorization criterion. 




Figure 4. Comparison of truncated models for p = 1. The continuous curve shows 
the case Z max = 1, i-e. the case of usual parallel update. The dotted curve shows 
( m = 2, given by (|18f) and the dashed curve shows the fundamental diagram for the 
shuffled update, obtained for Z max = oo. 



7. Application to pedestrian dynamics 

As mentioned in the beginning, the shuffled update has an interpretation ^21 EH! as 
an update procedure in two-dimensional cellular automata to avoid conflicts (situations 
in which more than one particle tries to access the same cell) occurring in parallel 
dynamics [T^j. An important example is the movement of pedestrians in a long 
corridor. To describe this movement more realistically we consider two straightforward 
generalizations of the model in the following. The simplest extension of the ASEP, 
which corresponds to a single-lane model is a model of several decoupled lanes. Here 
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pedestrians are not allowed to change lane, i.e. their y-coordinates are fixed. If in 
the initial state the pedestrians are distributed stochastically on sites (and thus lanes), 
their density in each lane is fixed and constant over time. The probability of finding i 
pedestrians in a particular lane is given by the hypergeometric distribution 




with mean p = N/ (LW) . The fundamental diagram of the model with decoupled lanes 

is |2Z| 

L 

J(p) = Hyp P w,L,L ( w-i)(0 • J(i/L), (20) 

i=0 

where J(p) is the single-lane fundamental diagram determined in Sec. El Using the 
approximation (|11|1. the result for a width of W = 10 cells and deterministic hopping 
(p = 1) in the thermodynamic limit is depicted in Figure El (a). One can see that the 
maximum is smoother and a little bit lowered in comparison to the corresponding Figure 
El (b). Note that lane- changing, an effect that is especially relevant for high densities, 
would lead to a reduction of the flow. 




Figure 5. Fundamental diagrams for two generalized processes with p = 1. The left 
diagram shows the result for the model with decoupled lanes for W — 10. The right one 
shows the case of increased velocity v max = 3 for W = 1 . The dashed curve represents 
the analytical result, the continuous curve the result from computer simulations. 

To reproduce the shift of the maximum flow to lower densities, observed in pedestrian 
dynamics, we introduce a larger maximum velocity f max , i.e. a larger number of cells 
that can be passed during one timestep [2H]. The simplest way is to allow a pedestrian to 
move the minimum of v max cells and the number of empty cells in front with probability 
p at each timestep. In the following we again consider the case W — 1 and p = 1 
only. For densities less or equal to l/(f max + 1) each of the pedestrians has at least t> max 
empty sites in front, i.e. the probability to find a particle in front is Po = and the 
flow is determinstically given as v max p. For higher densities it may happen that e.g. a 
pedestrian occupies a cell at time t+1 that has been occupied by a different pedestrian 
belonging to the cluster in front at time t. However we found in computer simulations 
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that this effect can be neglected (for relatively small v max ). Thus one can again describe 
the dynamics in the picture in which all clusters are updated in parallel. The system 
tries to separate the clusters by exactly t> max holes. The probability to have at least one 
hole in front is in this case (1 — p)/v max and the ratio of particles which do not have 
a particle in front is (p — (1 — p)/v Taax )/p. This defines the probability Pq, which can 
finally be written as [(f max + l)p — 1] 0{p — VC^max + l))/(^maxP)- With the use of (fT3*|) 
we obtain 

cP, for p < 1/Omax + 1), 



-'max/ 7 ) 



J(P, V max ,P = 1) = < V max p(l—p) I" (-max+l)p-l 1 (21) 



£ "maxP — 1 



else. 



1 + ^max)p - 1 

In Figure El (b) the resulting fundamental diagram in case of v max = 3 and p — 1 is 
depicted. The dashed line shows the analytical in comparison to the numerical result. 
In opposite to the case of maximum velocity 1, depicted in Figure El (b), the maximum 
of the fundamental diagram and the critical point (l/(v max + 1), pv m3iX ) do not coincide; 
this is qualitatively reproduced by the analytical result. However, to recover the exact 
result one would have to take longer ranged correlations into account. 



8. Discussion 



In this paper the asymmetric exclusion process (ASEP) with periodic boundary 
conditions and shuffled dynamics was studied. Using a mean-field approach, steady 
state properties like the fundamental diagram and distribution functions are derived 
approximatively that are in very good agreement with data from Monte-Carlo 
simulations of large systems. 

The model can be mapped onto a generalized zero-range process (GZRP) with usual 
parallel update, by interpreting the holes of the ASEP as the sites of the GZRP. These 
sites are occupied by the particles which follow the corresponding holes of the ASEP 
to the left. Further, the model can be mapped onto a zero-range process (ZRP) with 
shuffled update by interpreting the particles of the ASEP as sites of the ZRP. These 
sites are occupied by particles, that correspond to the holes of the ASEP. 

In the case p = 1 the exact results for the steady state distribution were presented. 
For densities p < 1/2 all particles are separated by at least one hole and thus the 
corresponding GZRP factorizes (see sec. EJ). Every configuration P(mi, . . . ,m,L-N) can 
thus be written as a product PfpQ~ 2N . For p > 1/2, no neighboring holes occur in 
the ASEP and thus in the corresponding ZRP sites are either occupied by exactly one 
particle or are empty with probabilities Pi = (l — p)/p and P = (2p — l)/p, respectively. 
However, since the shuffled dynamics is intrinsically stochastic, already for p — 1 the 
fundamental diagram is nontrivial and asymmetric. For p < 1/2 all particles can move 
deterministically and the stationary flow is identical to the deterministic limit of the 
parallel update which is known exactly. For higher densities p > 1/2 the probability 
for a particle to move depends on the number of particles in front and therefore the 
flow depends on the cluster-size distribution. However, mean-field theory yields a good 
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approximation for the flow. The two regimes are separated by a second order phase 
transition. In the case p < 1 the two regimes can no longer be distinguished and the 
flow is determined by the stochasticity for all densities. The fundamental diagrams 
become smooth, but still do not exhibit a particle-hole symmetry. The free flow and 
jammed regimes are no longer separated by a phase transition. It could be shown that 
in this more general case neither the GZRP, nor the ZRP have a factorized steady 
state. However, the mean-field theory yields again a very good approximation. In the 
limit p — > the model can be described by a product measure since the steady state 
distribution converges to the result of random-sequential updating. 

In contrast to the other updates investigated so far it is essential that we are dealing 
with distinguishable particles. The numbering is needed in order to define the update. 
This is not the case for random-sequential, parallel or ordered-sequential updates. Note 
that in the random-sequential case sites could be updated instead of particles without 
changing the stationary state. For the ordered-sequential dynamics only one particle 
has to be marked, namely the first one of the sequence. 

A truncated process in which at most two particles of a cluster can move and its 
exact solution for p = 1 was presented. It could be shown that this model factorizes 
into particle-cluster probabilities but not into hole-cluster probabilities. Since in shuffled 
dynamics the probability that the third particle in a cluster moves is relatively small, 
the truncated dynamics gives also a good approximation for the model with shuffled 
update. 

It has been shown ^2j that cellular automata with shuffled update have an 
interpretation as pedestrian flow models. In this sense the ASEP with shuffled update 
models the directed motion of pedestrians in a long corridor of small width. To 
allow pedestrians to walk side by side, we generalized the model to a two-dimensional 
scenario with decoupled lanes, in which pedestrians can not change lanes. Further 
a generalization to higher maximum velocities was given. For both generalizations 
approximations were presented that are in good agreement with Monte-Carlo data. 

We want to mention that the idea of ordered sequential updates (forward or 
backward ordered sequences) can be generalized 29J: Consider a general ordered 
sequential update with sequence n = tt(2), . . . , ir(N)), which now is quenched 

for all timesteps. 7r(l) denotes the update number of the particle with number 1 and so 
on, where the particles are numbered from left to right. The forward sequential dynamics 
is the special case of Tip = (1,2,..., N) and the backward sequential dynamics is the 
special case of ttb = (N, N — 1, . . . , 1). As far as we know, other update sequences 
have not been studied yet. However, it is easy to show that for N even all ((iV/2)!) 2 
sequences, in which the even (odd) particles are updated first, and then the odd (even) 
particles, have a factorized steady state distribution j21j. Note that these updates are 
equivalent to a sublattice-parallel-like update, in which in the first half of a timestep 
all even (odd), and in the second half, all odd (even) particles are updated in parallel. 
However, the usual ordered sequential update does not factorize with a general update 
sequence. Hence, it would be interesting to find a condition for a model with sequence 
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n (and transfer matrix hP^ . . . h^ N ') to have a factorized steady state. Among that, if 
an update with sequence n generates the steady state F(tt), it is interesting to consider 
the average of the steady states generated by all permutations of the permutation group 
S: (F) = F(tt)/N\. We found in computer simulations that this steady state 

differs from the steady state of the shuffled update distinctly, what can not be trivially 
anticipated. 

Summarizing, maybe the most surprising result is the fact that the ASEP with 
shuffled update does not factorize in contrast to the updates investigated before. Here 
the distinguishibility of the particles seems to be important whereas in the updates 
considered previously the particles are basically indistinguishable. This points warrants 
further investigations. 
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